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ABSTRACT 

We perform collisionless V-body simulations to investigate the evolution of the structural and kinematical 
properties of simulated thick disks induced by the growth of an embedded thin disk. The thick disks used in 
the present study originate from cosmologically-common 5:1 encounters between initially-thin primary disk 
galaxies and infalling satellites. The growing thin disks are modeled as static gravitational potentials and we 
explore a variety of growing-disk parameters that are likely to influence the response of thick disks. We find 
that the final thick-disk properties depend strongly on the total mass and radial scale-length of the growing thin 
disk, and much less sensitively on its growth timescale and vertical scale-height as well as the initial sense of 
thick-disk rotation. Overall, the growth of an embedded thin disk can cause a substantial contraction in both 
the radial and vertical direction, resulting in a significant decrease in the scale-lengths and scale-heights of 
thick disks. Kinematically, a growing thin disk can induce a notable increase in the mean rotation and velocity 
dispersions of thick-disk stars. We conclude that the reformation of a thin disk via gas accretion may play a 
significant role in setting the structure and kinematics of thick disks, and thus it is an important ingredient in 
models of thick-disk formation. 

Subject headings: cosmology: dark matter — cosmology: theory — galaxies: formation — galaxies: kinemat- 
ics and dynamics — galaxies: structure — methods: numerical 



1. INTRODUCTION 

Hierarchical models of cosmological structure formation, 
such as the c urrently favor ed cold dark matter ( "CDM ) 
paradigm (e.g., Whi te & Reeslfl978l: iBluinenthal et alj[l98l . 
generically predict that galaxies are built via the con- 
tinuous accretion of smaller systems. Indeed, a grow- 
ing body of observational evidence has recently con- 
firmed this prediction with the discovery of tidal streams 
and complex stellar structures in the Milky Way (MW) 
(e.g.. Ilbata et alj[l994t iHelmi et alj|1999t lYannv et al.1 2000j; 
Ibata et alj2001bHNewberg et al.ll2002l;lMaiewski e7ai]2003l: 
Martmez-Delgado et alj2005UBelokurov et al.l2006l). the An- 
dromeda galaxy (Ilbata et alj l2001a|F erguson et alj 120021 
l2005t iKahrai et al.ll2006t Ilbata et alj|2007b. and beyond the 
Local Group (e gjMalin & Hadleyl | 1997HShang et al| ll99l 
iPeng et alJEOOllForbes et al.ll2003t|Pohlen et alj|2004b . 

In the context of CDM, bombardment by infalling satellites 
may lead to the heating of thin galactic disks and the subse- 
quent f ormation of t hick d isk s (e.g.. lOuinn & Goodman! 
19861 IWalker et aU [l996l iVelazauez & White! 11999b 
lHavashi & Chibal l2006t IVillalobos & Helmil 120081 here- 
after lPaper It iKazantzidis et alj|2008t IRead et al.ll2008l). an d 
can even cause thin-disk destruction (Purcell et al. 2009). 
Depending on when the most substantial accretion event has 
occurred, a new thin disk may reform, for ex ample, via the 
coolin g of hot gas in t he gala ctic halo (e.g., IWhite & Reesl 
[19781: IMo et al.lll998l: iBaughl l2006h and be observable at 
present. 

In this scenario, the dynamical effects of the cooling gas 
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on the properties of the post-accretion thickened stellar distri- 
bution could be significant. Indeed, as the gas cools down 
and slowly accumulates at the center of the system, it can 
induce concomitant contraction of the heated stellar com- 
ponent due to its gravity (see e.g., lElmegreen & Elmegreen 
I2006t IKazantzidis et al.l 120091) . This phenomenon is similar 
in spirit to the contraction of dark matter halos during the 
process of bar yonic cooling that leads to the formation of 
galaxies (e.g.. IZeldovich et alj[l980t iBarnes & White! 11984 
Blumenth al et al.lll986l) ~ 

Significant theore tical effort, in cluding b oth semi- 
analy ti c modeling (iToth & Ostrikerl 119921: iBenson et al.l 
120041: iHopkins et alj |2008|) and numerica l simulations 
(lOuinn & Goodman! 119861: lOuinn et all 119931: IWalker et al l 
1 996t iHuang & Carlbergl Il997t ISellwood et all 119981: 
Velazq uez & Whitd 119991 iFontetalJ 120011: lArdietalJ 



20031: iGauthier et al 
Kazantzidis et al. 



12006T lHavashi & Chibal l2006t 

120081 IRead etal.1 120081 iPurcell et all 
2009t IKazantzidis et al.1 120091: iMoster et all 120091) has been 
devoted to exploring the dynamical heating of galactic disks 
via satellite accretion events. However, very few studies have 
investigated the response of the resulting thickened stellar 
component to the gas accumulatio n and subsequent refor- 
mation of a thin-d isk (see, however. IKazantzidis et alJl2009t 
IMoster et al.ll2009l) . A full model of satellite-disk interactions 
in a cosmological context including gas dynamics and star 
formation would be ideal to obtain a complete picture of this 
physical process. Unfortunately, due to the fact that they are 
highly complex and very costly in terms of computational 
power, cosmological simulations that could address in a 
systematic way the structural and kinematical evolution of 
thick disks i n response to growing t hin-disks do not yet exist. 

Recently, IKazantzidis et alj (120091) performed dissipation- 
less jV-body simulations to investigate the influence of a 
slowly growing thin disk on the properties of an initially- 
thick galactic disk. These authors considered three models 
of growing disks with masses equal to 10%, 50%, and 100% 
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of the mass of the thick disk, and showed that the latter con- 
tracted vertically as well as radially in response to the growth 
of the thin-disk component. Not unexpectedly, the magnitude 
of the thick-disk structural evolution was found to depend 
sensitively on the total m ass of the growing disk. However, 
iKazan tzidis et al. (2009) focused on a fairly small region of 
parameter space and, in addition, they utilized a fully-formed 
thick-disk galaxy instead of a galactic model whose thick disk 
formed self-consistently by encounters between an initially- 
thin disk and infalling satellites. 

In the present paper, we expand upon this initiative by per- 
forming a systematic numerical study aiming to elucidate the 
effects of growing thin disks on both the morphology and 
kinematics of realistic thick disk s. For th e latter, we adopt a 
subset of the models presented in lPaper j . Given the complex 
interplay of effects (e.g., gas cooling, star formation) relevant 
to the formation and evolution of spiral galaxies and the out- 
standing issues regarding disk galax y formation in CDM cos- 
mogonies (e.g., Maver et al. 2008), we restrict ourselves to 
modeling the gravitational potential of a thin disk that slowly 
grows in mass over time using an ensemble of collisionless 
iV-body simulations. 

Our simulation set is carefully designed to permit an in- 
vestigation of a large parameter space and we explore several 
aspects of the growth of a thin disk that are likely to have 
an effect on the morphological and kinematical evolution of 
thick disks. In particular, we vary the growth timescale, final 
mass, scale-length, scale-height, and orientation of the grow- 
ing disk total angular momentum vector. In addition, we per- 
form experiments where the same thin disk grows inside thick 
disks that rotate in either a prograde or retrograde sense with 
respect to their halos. The present work establishes that the 
reformation of a thin disk may play a significant role in set- 
ting the structural and kinematic properties of thick-disk stars, 
and thus it is an important ingredient in models of thick-disk 
formation. 

The outline of this paper is as follows. In Section [2] we 
briefly describe the models adopted and the setup of the nu- 
merical experiments performed in the present study. Section[3] 
contains the results regarding the thick-disk morphological 
and kinematical evolution induced by the growth of a thin- 
disk component. A detailed comparison between the prop- 
erties of the initial and post-growth thick disks is presented 
in Section 2] Lastly, Section [5] summarizes our main conclu- 
sions. Throughout this work we use the terms "thick disk" 
and "heated disk" interchangeably to indicate the perturbed 
stellar component that resulted from the original interaction 
between a thin primary disk galaxy and an infalling satellite, 
and consists of stars from both the stellar disk of the primary 
galaxy and the disrupted system. 

2. METHODS 

We begin this section by briefly descr ibing the primary disk 
models used in this study. In Section l2~2l we introduce the 
thin-disk models and discuss the method we empl oy to grow 
them inside the thick disks. Lastly, in Section l2~3l we present 
the numerical experiments with growing disks of the adopted 
simulation campaign. 

2. 1 . Models of Thick Disks 

In lPaper fl, 25 dissipationless A^-body simulations were car- 
ried out to study the general properties of thick disks formed 
by single encounters between a primary thin disk galaxy and a 
massive satellite system. These simulations explored: (1) two 



models for the primary disk galaxy, whose structure and kine- 
matics resembled those of the MW at present ("z=0" experi- 
ments) and at z=l ("z=l" experiments); (2) two morphologies 
for the stellar component of the satellite, spherical and disky; 
(3) two mass ratios between the satellite and the primary disk 
galaxy (10% and 20%); and (4) three initial satellite orbital 
inclinations (i = 0°,30°, and 60°), in both prograde and ret- 
rograde directions with respect to the rotation of the primary 
disk. 

Both the primary disk galaxy and the satellite were ini- 
tialized as "live" A^-body models consisting of a dark mat- 
ter and a stellar c omponent. The dark halos of both sys- 
tems followed the iNavarro et alj J1997I) density profile and 
have been adiabatically contracted to acco unt for the pres- 
ence of the stellar component (see e.g., [Blumenthal et al. 
1986; Mo et al. 1998). The stellar component of the spherical 
satellite was set up with structure and kinematics consistent 
with the observed fu ndamental plane of dE+dSphs galaxies 
dde Rijcke et al.ll2005l) . On the other hand, the disky satellites 
have been constructed in the same way as the primary disk 
galaxies. Their scale-lengths and scale-heights are smaller 
than those of the primary disks, following the proportion be- 
twee n the v irial radii of the corresponding dark matter halos 
(see iPaper II for details). The orbital parameters of the en- 
counters were drawn fro m studies of in falling s ubstructures in 
a cosmological context dBensonll2005l see also iTormeiill 1 9971 ; 
iKhochfar & Burkerti 2006), and the satellites were released far 
away from the centers of the disk galaxies, ~ 35 (~ 50) times 
the disk scale-length for "z=0" ("z=l") experiments. 

The simulations of thick-disk formation were evolved for 
5 Gyr (4 Gyr) in the "z = 0" ("z = 1") experiments and the 
satellites have typically merged by t = 3 Gyr (f = 2 Gyr). 
By the final time, the systems are relaxed. The simulations 
were carried out with 5 x 10 5 particles in the host dark matter 
halo, and 10 5 particles in each of all the other components. 
The softenings of thic k-disk st ars range from 12 pc to 70 pc 
(we refer the reader to IPaper II for a complete description of 
the experiments), and the softening of growing disk particles 
corr espo nd to 1/5 of the growing disk scale-height (see Sec- 
tion |23). 

The primary aim of the present work is to study the most 
generic features of the evolution of thick disks in response 
to the growth of an embedded thin disk. For this reason, 
we sele ct a subset of t he thick-disk models presented in 
IPaper II (see Section |231 below). In particular, we focus on 
models that were produced via an encounter with a satel- 
lite on a low/intermediate orbital inclination. Assuming a 
dynamical heating mod el for the origin of thick disks (e.g. 
IKazantzidis et al.ll2008l) . such accretion eve nts are re levant to 
the formation of the thick disk of the MW (Paper I). This is 
based on the possible existence of a vertical gradient in the ro- 
tational velocity of the Galactic thick disk and on t he observed 
value of q- 7 Ior in the sola r neighborhood (e.g. iGirard et ail 
120061; IVallenari et"aill2006l) . In addition, due to the fact that 
satellite structure is not es sential in establishing the structural 
properties of thick disks (Paper I), we only consider thick- 
disk models that were produced by interactions with spherical 
satellites. 

2.2. Models of Growing Disks 

The growing disks follow an exponential distribution in 
cylindrical radius R, and their vertical structure is mod eled in 
the standard way as a collection of isothermal sheets ( Spitzer 
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TABLE 1 

Summary of Numerical Experiments 



Growing thin disk Thick disk 







Rd 


Zd 


T 








Run 


M d /M thw k 


(kpc) 


(PC) 


(Gyr) 


alignment 


rotation 


inclination 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


Reference 


5 


3 


125 


1 


halo 


prograde 


30° 


model 














Al 


5 


3 


125 


1 


halo 


prograde 


0° 


A2 


5 


3 


125 


1 


halo 


prograde 


60° 


Bl 


5 


3 


125 





halo 


prograde 


30° 


B2 


5 


3 


125 


5 


halo 


prograde 


30° 


CI 


2 


3 


125 


1 


halo 


prograde 


30° 


C2 


5 


3 


25 


1 


halo 


prograde 


30° 


C3 


5 


1 


125 


1 


halo 


prograde 


30° 


D 


5 


3 


125 


1 


thick disk 


prograde 


30° 


El 


5 


3 


125 


1 


halo 


prograde 


0° 


E2 


5 


3 


125 


1 


halo 


retrograde 


0° 



NOTES. — Column (1): Labels for the numerical simulations. Column (2): Mass of the growing thin disk in units of the initial thick disk. Column (3): Radial 
scale-length of the growing thin disk in kpc. Column (4): Vertical scale-height of the growing thin disk in pc. Column (5): Growth timescale of the thin disk in 
Gyr. Column (6): Alignment of the rotation axis of the growing thin disk. Column (7): Sense of rotation of the thick disk with respect to its halo. Column (8): 
Initial orbital inclination of the infalling satellite that produced the initial thick disk. 




FIG. 1 . — Initial (upper panels) and final (bottom panels) surface brightness contours of thick-disk stars viewed face-on (left panels) and edge-on (right panels) 
in the simulations with growing thin disks. Results are presented for experiment A with a thin-disk growth timescale of r = 1 Gyr and the orbital inclination of 
the accretion event that formed the initial heated disk is indicated in each panel. Contours correspond to equal surface brightness levels (0.39 mag arcsec ) in 
the V-band for face-on (20.25-25.35 mag arcsec -2 ) and edge-on (19.09-25.35 mag arcsec -2 ) views. Thick disks experience strong radial and vertical contraction 
in response to the growth of a new thin disk component. 



1942): 

P Mn(R,z) = Kexp (-^j sech 2 (j^j , (1) 

where M d , Rj, Zd denote the mass, radial scale-length, and 
(exponential) vertical scale-height of the disk, respectively, 
and K = M d / SirRdZd ■ 

Each growing-disk simulation was performed using the fol- 
lowing procedure: (1) insert a massless Monte Carlo particle 
realization of the desired disk model inside the primary disk 
galaxies. (2) increase the mass of this distribution to its final 
value linearly over a timescale r, according to the following 
law: 

M{t) = M d {t/r) 0<t<r. (2) 



During the growth period, the growing disk remains rigid 
and their particles are fixed in place (i.e., behaving like a fixed 
disk potential), while the "live" particles composing the thick 
disk and its dark matter halo are allowed to achieve equi- 
librium as the mass of the system grows. Throughout the 
experiments, all other properties of the growing disk (e.g., 
scale-length, scale-height) are kept constant. All numeri- 
cal simulations of growing disks were carried out w ith the 
multi -stepping, parallel, tree N-body code PKDGRAV (IStadell 
l200l . 

2.3. Numerical Experiments with Growing Disks 

Overall, we have performed 1 1 simulations of the growth 
of a new thin disk within a subset of the thick disks presented 
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in iPaper H We have explored a variety of parameters related 
to the growing disks that may affect the thick-disk evolution, 
including: (1) the growth timescale: r = Gyr, 1 Gyr, and 
5 Gyr 4 ; (2) the final mass: Md = 2M t hi c k and 5M t h; c k, where 
denotes the mass of the initial thick disk 5 ; (3) the verti- 
cal scale-height: Zd = 25 pc and 125 pc 6 ; (4) the radial scale- 
length: Rd=l kpc and 3 kpc; (5) the orientation of the angular 
momentum vector: aligned either to the angular momentum 
of the dark matter halo or to the angular momentum vector of 
the thick disk; and (6) the sense of rotation of the initial thick 
disk with respect to its halo: prograde or retrograde. Lastly, 
we also compare the evolution of initial thick disks produced 
by infalling satellites on different orbital inclinations, i = 0°, 
30°, and 60°. 

The runs are labeled as A, B, C, D, and E and summarized 
in Table Q] Our reference model is characterized by the fol- 
lowing set of parameters for the growing thin disk: r = 1 Gyr, 
Md = 5M t hi c k, Zd = 125 pc, Rd = 3 kpc, and a thick disk that 
is the result of a (prograde) satellite accretion event with an 
initial orbital inclination of i = 30°. In what follows, we refer 
to experiment "X" to denote comparisons between Run "X" 
and the reference model, except for run E. 

All numerical experiments of growing disks correspond to 
a total time of t = t + 0.5 Gyr (except for the case of instanta- 
neous growth which is evolved for 1 .5 Gyr to match the simu- 
lations with a growth timescale of r = 1 Gyr). Apart from the 
production simulations described above, we performed two 
additional sets of experiments. First, we evolved a subset of 
the initial thick disks in isolation for 1 .5 Gyr and confirmed 
that their properties remained fairly unmodified. Second, after 
the full growth of the new thin disks, we evolved the resulting 
composite disk galaxies for another 4 Gyr in isolation, con- 
firming that their properties do not evolve appreciably with 
time. Therefore, we conclude that the evolution of the thick- 
disks structural and kinematic properties presented next is due 
solely to the adiabatic compression caused by the growth of 
the thin disks. 

3. RESULTS 

3.1. Structural Evolution of Thick Disks 

3.1.1. General Features 

In general, the growth of a new thin disk induces three char- 
acteristic structural changes in all initial thick disks. First, it 
causes significant contraction in both the radial and vertical 
direction. Figure Q] focuses on experiment A and shows ini- 
tial and final surface brightness contours of thick-disk stars 
viewed both face-on and edge-on. We stress that no stars from 
the growing thin disks are included in this analysis. The fact 
that the final contours appear much closer together and signif- 
icantly less extended in the outer regions illustrates the con- 

4 Note that these values bracket the range of potential timescales for the 
disk growth including strongly non-adiabatic and adiabatic ones. While 
adopting an instantaneous growth for the disk is obviously unphysical, it still 
represents an interesting limiting case for the disk growth timescale. 

5 Note that the value of Md = 2M,hi c k = 2.8 X 10 10 Mq is within the range 
of the total gas mass accreted by the MW disk since z ~ l (1.5— 3 X 1() 10 Mq 
assuming a constant infall rate of 2— 4 Mq / yr) according to an alytical models 
for the evolution of the MW disk in a cosmological context INaab & Ostrike3 
2006). Although there are still large uncertainties in this estimate, the larger 
value of Mj = 5M t i,i c k ls relevant to a syste m whose thin-thic k mass ratio is 
similar to that of the Galaxy at present (e.g., June et al. 2008) 

6 Note that these values are consistent with the scale-heights of 
known, young, star-forming disks observed b oth in external galax- 
ies (e.g.. |Wainscoat et al . 1989; Matthews 2000) and in the MW (e.g., 
BahcalT & Soneira 1980; Reid & Majewski 1993). 



traction experienced by the disk stars. At the same brightness 
level, the innermost contours of the final systems are more 
extended compared to those of the initial ones. This indi- 
cates that the surface brightness associated with the central 
regions of the final systems increases. It is also important 
to note that the face-on surface brightness contours become 
rounder after the growth of the thin disk. This confirms re- 
sults of earlier studies regarding the effects of bar yonic dis- 
sipation on the shapes of da rk matter halos (e.g., [Dubinski 
[l994llKazantzidis et al.ll2004l) . 

Figure [2] compares for the same experiment the initial ra- 
dial positions of initial heated-disk stars to those after the 
growth of the new thin disk. Only stars within R < 15 kpc 
and \z\ < 1 kpc in the final contracted thick disks are included 
in this analysis. The figure shows that the relation between 
initial and final radius is remarkably linear in all cases with 
practically the same slopes, regardless of the initial orbital in- 
clination of the infalling satellite. In addition, the structural 
contraction induces a significant migration of thick-disks stars 
inwards. Indeed, stars that were initially located at R ~ 14 kpc 
end up at R ~ 8 kpc after the growth of the thin disk. 

The angular momenta of individual thick disk particles are 
well conserved after the growth of the new thin disk, as it 
is shown in Figure [3] for experiment Al. In general this is 
found in all experiments, except in Bl, as it is expected given 
its non-adiabatic nature. Particles within R < 3 kpc present 
a less clear conservation of their angular momenta, which is 
probably due to the exchange of angular momentum between 
halo particles and thick disk particles in the inner region. 

The second common feature seen in all experiments is that 
the col d components present in the initial heated disks (see 
IPaper J) are preserved after the thin-disk growth. As we dis- 
cuss in the next section, the mass fraction associated with 
these components remains below ^25% of the total mass of 
the heated disks in all experiments. In principle, these thin 
components could represent old thin-disk populations given 
their similar spatial distribution and kinematics. 

As discussed in Paper I, the satellite accretion events which 
produce the thick disks create flares in the outer region s of the 
primary disk galaxies (see also Kazantzidis et al. 2008). How- 
ever, in spite of such flares, it was also found that (radially 
integrated) vertical surface brightness profiles of thick disks 
are very well represented by two sech 2 components. Another 
common characteristic to all experiments described here is the 
fact that the growth of a new thin disk also changes the vertical 
structure of thick disks, especially in the outer regions. The 
relevant analysis is presented in Figure H] which shows results 
for run Al. This figure compares the final vertical surface 
brightness profile (integrated for R < 10 kpc) of the contracted 
thick disk to that of the initial thick disk evolved in isolation 
(for the same amount of time), and demonstrates that they 
are fairly different. More specifically, the growth of the thin 
disk causes the vertical surface brightness profiles to become 
substantially narrower (a closer inspection shows a significant 
contribution of satellite stars in the region \z\ > 1.5 kpc, es- 
pecially at larger radii). While a two-component sech 2 de- 
composition provides an accurate description of the vertical 
structure of the initial thick disk down to fairly low surface 
brightnesses, the growth of the new thin disk highlights the 
need for a more complicated functional form to describe the 
final vertical disk structure at large heights. 

3.1.2. Experiments from A to E 
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FIG. 2. — Initial and final cylindrical radial distance of heated disk stars before and after the structural contraction. Results are presented for experiment A. The 
dashed line shows the best linear fit. For clarity, only heated disk stars are shown. 




FIG. 3. — Ratio of final-to-initial angular momenta of particles belonging to 
the isolated thick disk (i.e., evolved without the growing new thin disk) (left 
panel) and the corresponding contracted thick disk in experiment Al (right 
panel). In both cases the thick disks have been evolved for the same amount 
of time. The dashed line indicates where |iftnal /-^-initial I = !• 
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FIG. 4. — Final vertical surface brightness profiles of the isolated thick disk 
(i.e., evolved without the growing new thin disk) (left panel) and the corre- 
sponding contracted thick disk in run Al (right panel). The former has been 
evolved in isolation for the same amount of time as the latter and profiles 
are obtained via integrating within R < 10 kpc. Dashed and solid lines show 
results for the best-fit with one- and two-component sech 2 decompositions, 
respectively. The growth of a new thin disk induces changes to the vertical 
structure of the initial thick disks, highlighting the need for a more compli- 
cated functional form to describe the final disk structure at large heights. 



Figures [5] to QT| present results for all experiments from A 
to E. In particular, these figures show (1) the vertical surface 
brightness profiles of the thick disks at various projected radii 
(including the best two-components sech 2 fits) before and af- 
ter the growth of the new thin disk; (2) the surface density de- 



composition (into cold/thinner and hot/thicker components) 
of the final thick disks, including the parts of the profiles con- 
sidered to estimate the respective scale-lengths and the mass 
fraction associated with the cold component; and (3) a com- 
parison between scale-lengths and scale-heights of thick disks 
(for cold and hot components) before and after the growth of 
the new thin disk. All of these properties are computed within 
R < 10 kpc and \z\ < 4 kpc. Lastly, unless otherwise explicitly 
stated, we do not include stars from the growing thin disks in 
the analysis presented below. 

Briefly, the decomposition of the vertical surface bright- 
ness profiles is obtained by fitting, at equally spaced pro- 
jected radii on the (properly aligned) thick disks, the func- 
tion L(z) = L\ sech 2 (z/2zi)+L2 sech 2 (z/2z2), where L, and Zi 
are the local brightness on the plane and local scale-heights, 
respectively, and both are free parameters (although in gen- 
eral zi < Z2). In these fits, stars from both the disk and the 
(disrupted) satellite are included within the spatial range men- 
tioned above. For each thick disk one global scale-height per 
component is finally obtained by averaging the local scale- 
heights over all radii, weighting by the local brightness on the 
plane. On the other hand, the decomposition of the surface 
density profiles of thick disks is obtained by simply consider- 
ing the region < \z\ < 0-5z\ as dominated by the cold/thinner 
component, and the region z\ < \z\ < 4 kpc as dominated by 
the hot/thicker component. Finally, the scale-lengths of each 
component are computed by applying a linear fit to In £(/?), 
avoiding non-axisymmetries associated to both the central re- 
gions an d the ve ry outskirts. For more details, we refer the 
reader to Paper II 

Experiment A (Figure|3J shows that the final mass fraction 
associated to the cold component is larger than the initial one 
for thick disks formed after an either low or intermediate incli- 
nation merger (i = 0° - 30°), as opposite to the case when the 
thick disk is formed after a high inclination merger (i = 60°), 
where the final mass fraction is smaller than the initial one. It 
is also found that in an absolute sense the final mass fraction 
associated to the cold component remains within the range 
~ 14-20%, as before the growth of the new disk. The excep- 
tion is the i = 60° case in which the mass of the cold remnant 
may amount up to 50% of that of the initial cold component. 
This fraction however is relatively poorly determined because 
it depends strongly on the region selected to fit the exponential 
profile to the surface density which appears to show a change 
of slope around R ~ 4 kpc. 

Experiment B (Figure© shows that the thick-disk structural 
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FIG. 5. — Experiment A. Final thick-disk properties induced by the growth of the fiducial thin disk on a timescale of r = 1 Gyr. The initial thick disks were 
formed by satellite accretion events with initial orbital inclinations of i = 0° (left panels), ;' = 30° (middle panels), and = 60° (right panels). From top to bottom: 
Rows 1 and 2 show the initial and final vertical surface brightness profiles of the thick disks. Row 3 shows results for final thick-disk surface density profiles 
decomposed into a "cold" and a "hot" component. The mass fractions associated with these components are indicated in each panel. Row 4 presents initial (filled 
symbols) and final (open symbols) scale-lengths and scale-heights for the "cold" and "hot" components as a function of the orbital inclination of the encounter 
that produced the initial thick disks. 
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evolution does not have a strong dependence on the growth 
timescale of the new thin disk for r=l and 5 Gyr. In both 
cases, the peak of the final vertical surface brightness profiles 
becomes narrower and is characterized by a central value that 
is ~ 1 magnitude brighter than that of the initial configuration. 
However, at distances R > 4 kpc, the surface brightness levels 
near the plane (z = kpc) are lower compared to the initial 
thick disk. Both features reflect the radial and vertical con- 
traction discussed above. Moreover, the initial scale-heights 
of the cold and hot components decrease by 20-35% after 
the accretion event. While the scale-lengths have decreased 
by ~ 30% for both the cold and hot component, the radial 
surface density profiles are very similar for both timescales. 

On the other hand, the most dramatic feature of the exper- 
iment with the instantaneous growth is the "break" at R ~ 
3-4 kpc of the originally exponential surface density profiles 
of both the cold remnant and the hot/thicker component. For 
t = Gyr the central surface densities reach similar values to 
those of the t = 1 Gyr and r = 5 Gyr cases, while in the region 
R > 4 kpc they show little evolution compared to the initial 
system. 



Experiment CI shows that a factor of 2.5 more massive 
growing thin disk causes roughly twice more contraction in 
the initial thick disk in terms of scale-length and scale-height. 
Interestingly, the mass fraction associated with the cold com- 
ponent of the thick disk increases by only ^3% (this equiva- 
lent to a mass increase in the cold component of roughly 10%) 
(Figure [7]). On the other hand, the thickness of the growing 
thin disk (experiment C2) does not seem to impact signifi- 
cantly the structural evolution of the thick disks, as shown 
in Figure [8] where the left and right panels correspond to 
Zd = 25 pc and Zd = 125 pc, respectively. For example, in 
both cases the scale-lengths decreased by ~ 30%, while the 
scale-heights decreased by ~ 40%. Finally, variations in the 
scale-length of the new disk (experiment C3) lead to a change 
in the mass fraction associated with the cold component of 
the thick disk. Thin disks with a scale-length 3 times smaller 
(Figure |9]l reduce the mass of the cold remnant to one third. 
This can be expected since a growing disk with a larger radial 
extent would be able to attract more mass onto the midplane 
across the system. No significant differences are detected in 
the final scale-lengths and scale-heights of both components 
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FIG. 8. — Experiment C2. Final thick-disk structural properties for different vertical scale-heights of the growing thin disk: zd = 25 pc (left panels) and 
Zd = 125 pc (right panels). Rows 1 to 3 as in Figure|5] 



in the contracted thick disks. However, the estimation of the 
final scale-length of the cold remnant in the case of a new disk 
with R t i = 1 kpc should be taken with caution, since it could 
be affected by the "break" feature of the surface density and 
by the region chosen to derive the scale-length. 

We remind the reader that the previous analysis does not 
include stars from the growing thin disks. Therefore, it is 
important to clarify the extend to which the structural de- 
composition of the final thick disks depends on the inclusion 
of the growing thin disks. Such an attempt is presented in 
Figure [9] which compares the final thick-disk vertical surface 
brightness profiles in experiment C3 (Rj = 1 kpc) with and 
without the inclusion of the growing thin disk. Not unex- 
pectedly, when the growing thin disk is included in the de- 
composition, the peak of the surface brightness profiles in- 
creases at all radii. However, the scale-heights of both "cold" 
and "hot"components are not substantially modified. Specif- 
ically, when the growing thin disk is included, scale-heights 
change from 140 pc and 585 pc to 130 pc and 650 pc. We 
find a similar behavior in the rest of our experiments for the 
scale-heights. It is important to note that other properties can 



be affected more by the inclusion of the growing thin disk in 
the analysis (e.g., scale-lengths). However, the magnitude of 
these differences depends sensitively on the criteria that are 
used to define the regions dominated by the "cold" and "hot" 
component, that is |z| < 0.5z co id and z co i d < |z| < 3 kpc, re- 
spectivel y, in term s of the scale-height of the "cold" compo- 
nent (see PaperJ). 

In the numerical simulations discussed in Paper I, the an- 
gular momenta of the dark matter halos and thick disks were 
often slightly misaligned by few degrees 7 . In experiment D, 
we utilized the thick disk produced by an accretion event with 
an initial inclination of i = 30°. The angular momentum vector 
of the adopted thick disk is misaligned by ~ 6° with respect 
to the halo angular momentum. For this experiment we have 
allowed the angular momentum of the growing thin disk to 
be aligned with the angular momentum of either the halo or 

7 Cosmological simulations show that the angular momentum of the bary- 
onic component of a galaxy correlates well with the angular momentum 
of the parent dark matter halos, with a typical misalignment of ~ 20° 
ISharma & S teinmetz 2003). 
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FIG. 9. — Experiment C3. Final thick-disk structural properties for different radial scale-lengths of the growing thin disk: Rj = 1 kpc (left panels) and Rj = 3 kpc 
(right panels). Rows 1 to 3 as in Figure|5] The top right panel shows the results when the growing thin disk is included in the decomposition analysis. 



the thick disk. Figure[lO]shows that no significant differences 
can be detected in the structure of the thick disks for either 
alignment of the growing disk, which is not surprising given 
the rather small misalignment between the initial thick disk 
and the halo angular momenta. In principle, the differences in 
the evolution of the properties of the thick disk could become 
larger for a more significant initial misalignment. 

In experiment E (Figure ITTb the adopted initial thick disks 
have been produced by an accretion event with an initial in- 
clination of i = 0°. We considered two cases where the in- 
falling satellite was on a prograde and on a retrograde orbit 
with respect to the rotation of the pre-existing galactic disk. 
In both cases, a significant amount of angular momentum 
is transferred from the satellite to the dark matte r halo (see 
IVitvitska et"alll2002t iHetznecker & Burkertll2006l) and much 
less into the pre-existing disk. The latter is because by the 
time the satellite has reached the disk the orbit is much smaller 
than the initial one, and it has lost most of its bound mass. As 
a result, the halo and the resulting thick disk rotate in the same 
sense after the accretion event in the prograde case, while in 
the retrograde case they rotate in opposite directions. We note 



that our initial dark matter halos are non-rotating, but the same 
conclusions would hold for halos with typical amounts of ini- 
tial rotation (see below). 

The surface density profiles in Figure QT| show that the 
growth of a thin disk on thick disks with opposite sense of 
of rotation induces essentially the same increase in the mass 
fraction associated with the thin remnant (from ^15% ini- 
tially to ~ 22%), and that the decrease in the scale-length of 
the thin remnant is of similar amplitude (~ 30%). The most 
significant difference seems to be in the scale-lengths of the 
thick disk, where in the prograde case it is induced a decrease 
of 25%, while in the retrograde one it is induced only a 7% 
decrease (as for experiment C3 this has some dependency on 
the region considered to estimate the scale-lengths). On the 
other hand, the evolution of the scale-heights of both the thin 
remnant and the thick disk do not show a significant differ- 
ence between prograde and retrograde cases. 

This experiment was motivated by recent observations of 
counter-rota ting thick disks (with respect to their thin disks) 
detected by Yoachim & Dalcanton (2005). It has been sug- 
gested that such counter-rotation may pose difficulties for the 
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— Experiment D. Final thick-disk structural properties for different alignments of the total angular momentum vector of the growing thin disk: either 
angular momentum of the dark matter halo (left panels) or along the angular momentum vector of the thick disk (right panels). Rows 1 to 3 as in 



disk heating scenario as a viable model for the formation of 
thick disks. However, an encounter between a pre-existing 
thin disk galaxy and a massive satellite on a retrograde orbit 
can, in principle, produce a thick disk and a counter-rotating 
dark matter halo. In order to illustrate this we can estimate 
the total angular momentum of a halo after such encounter as 
L]? al ° = Lf dlo + Ll a \ where the last term on the right is the initial 
orbital angular momentum of the infalling satellite. Lj 13 ' can 
be computed using the halo spin parameter, A, and its virial 
properties as Lf dl ° = %/2AM vil -V vir fl vir dBullocket al.l2001l) . On 
the other hand, L? at = M^RqVq, in terms of the satellite mass 
and its initial orbital radius and velocity. In the case of ex- 
periment E2, M vir = 5.07 x 10 u Mq, V vir = 133.87 km s" 1 , 
R vir = 122.22 kpc, M sat = 10 n M Q , R Q = 83.9 kpc, and 
Vb = -78.8 km s -1 (negative sign indicates a retrograde or- 
bit). Thus, assuming a typical value for the spin parameter, 
A = 0.035, we have Lf l ° = -2.5 x 10 I4 M Q kpc km s" 1 , where 
the negative sign shows that the halo is counter-rotating with 
respect to its initial state. 



Freshly accreted cold gas from the galactic halo would have 
its angular momentum aligned with that of the dark halo and 
therefore could form a young thin disk that is counter-rotating 
with respect to the older thick disk. In this case, the new disk 
would also counter-rotate with respect to any old thin remnant 
present. 

3.2. Kinematical Evolution of Thick Disks 
3.2.1. General Features 

Similarly to the structural evolution, the growth of a thin 
disk induces two characteristic changes in thick-disk kine- 
matics that are common to all experiments. In particular, a 
growing thin disk causes an increase in both the mean rota- 
tional velocity and the velocity dispersions of the thick disks. 
Both effects can be explained by the fact that the growing thin 
disk adds mass to the system and by the subsequent structural 
contraction described above. 

3.2.2. Experiments from A to E 
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FIG. 1 1 . — Experiment E. Final thick-disk structural properties for different sense of rotation of the initial thick disk: prograde (left panels) and retrograde 
(right panels). Rows 1 to 3 as in Figure|5] 



Figure Q~2] presents the initial and final global kinemat- 
ics of the simulated thick disks as a function of galactocen- 
tric radius. We characterize the thick-disk kinematical re- 
sponse to the growth of a new thin disk through the evolution 
of the mean rotational velocity, (V^), and velocity ellipsoid 
(ojfjCr^jCTj), where <tr, cr^, and cr z correspond to the radial, az- 
imuthal, and vertical velocity dispersions, respectively. These 
properties are computed within R < 10 kpc and \z\ < 4 kpc 
and as before no stars from the growing thin disks have been 



included in the analysis. It is important to mention that we do 
not find significant differences between the kinematics of the 
global thick disks (i.e. "cold" + "hot" components) and those 
from the "hot" component alone in the case of <jr and o^. 
However, ov and of the "hot" component are larger, as ex- 
pected, by ~ 20 kms -1 and ~ 40 kms" 1 , respectively (mainly 
at smaller radii). 

The results of experiment A demonstrate that after the 
growth of the new thin disk, the velocity dispersions of the 
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final thick disks maintain the original trends with the satellite 
orbital inclination (that is, the hotter the initial thick disk the 
hotter the final one). We note that the mean rotational veloc- 
ity in the 60° case shows a clear decrease at larger radii which 
is due to the increasing skewness of the total V<j, distribution 
towards lower velocities. This is a consequence of the contri- 
bution of the accreted stars, which are preferentially found in 
the outskirts and which rotate more slowly with respect to the 
stars of the pre-existing thin disk. 

The findings of experiment B indicate that for r=l and 
5 Gyr the global kinematics of the final thick disks do not de- 
pend on the growth timescale of the new disk. Interestingly, 
the effect of the instantaneous growth is more complex. In this 
case, the mean rotation of the thick disk is not substantially af- 
fected by the growth of the new disk (except by some decrease 
in the inner region). This is consistent with the corresponding 
surface density profile (Figure|6]l which shows little variation 
with respect to the initial configuration (except by a distinc- 



tive "break" at R ~ 4 kpc). Although a significant amount of 
mass has been added to the system via the growth of the new 
thin disk, this growth was strongly non-adiabatic. This fact 
implies that the eccentricity of the stellar orbits in the thick 
disk must have changed significantly, or equivalently, that the 
amount of ordered (circular) motion must be much smaller. 
Thus in practice, the extra potential energy due to the new thin 
disk must have gone into random motions in the plane as can 
be seen from the central panels in Figure Q~2] which clearly 
show an increase in the velocity dispersions, particularly at 
R ~ 4 kpc. 

The results of experiment CI show that a more massive 
growing disk induces a larger increase in both the mean ro- 
tation and the velocity dispersions of the final thick disks. By 
resorting to the Jeans equations (Binnev & Tremaine 119871 
Eq. 4.29), it may be expected that a thin disk roughly twice 
as massive causes an increase of a factor of 2 in all thick-disk 
velocity dispersions and in the mean rotational velocity at a 
given radius. 

The findings of experiment C2 suggest that for our choices 
of vertical scale-heights, id = 25 and 125 pc, the kinematical 
evolution of thick-disk stars does not depend essentially on 
the thickness of the growing disk. Of course, this conclusion 
is valid as long as the mass of the growing disk is the same as 
in our case. On the other hand, the results of experiment C3 
highlight that at a fixed mass a more radially compact growing 
disk is responsible for generating a larger increase in both the 
mean rotational velocity and velocity ellipsoid in the inner 
regions of thick disks. 

The results of experiment D show that the relative align- 
ment of the new disk (with respect to either the halo or the 
initial thick disk angular momentum) does not seem to affect 
the kinematical evolution of thick disks. This could also have 
been expected from the lack of structural evolution shown in 
Figure [TO] Note, however, that the misalignment between the 
halo and initial thick disk angular momenta is only ~ 6°. Such 
misalignment may be possibly too small to induce significant 
differences in both the structure and kinematics of the final 
thick disks. 

The findings of experiment E demonstrate that the evolution 
of all three velocity dispersions is very similar for both pro- 
grade and retrograde thick disks, as expected from Figure fTTI 
Note that for the retrograde case there is a significant decrease 
in the mean rotation of the thick disk at R > 4 kpc. Simi- 
larly to experiment A2, a closer inspection shows that there 
is significant skewness in the distribution of towards lower 
velocities. This is due to the contribution of satellite stars at 
large radii which rotate more slowly than disk stars. Overall, 
in the retrograde case the evolution in the kinematics of the 
thick disk is found to be very similar to that in the prograde 
case. The most significant difference is found in a z where the 
retrograde case leads to a somewhat higher dispersion at large 
radii. 

4. COMPARISON WITH INITIAL THICK DISKS 
4.1. Structure 

The results presented in Figure 8 of iPaper 3 show that the 
outer isophotes of the simulated thick disks measured at sur- 
face brightness [iy > 25 mag arcsec -2 , namely > 6 mag be- 
low the central peak /iy ~ 19 mag arcsec" 2 , are consistently 
more boxy than the inner ones. This characteristic feature is 
present in all the thick disks simulated in IPaper 1 implying 
that the detection of such a degree of boxiness in real galaxies 
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FIG. 14. — Final mean azimuthal velocity of thick-disk stars as a function of galactocentric distance, R. Results are presented for experiment A at various 
scale-heights above the disk plane. The vertical gradient of the mean azimuthal velocity depends sensitively on the orbital inclination of the infalling satellite that 
produced the initial thick disks. 
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FIG. 15. — Final velocity dispersion ratio ct./ctr of thick-disk stars as a 
function of galactocentric distance, R. Results are shown for experiment A. 
Even after the growth of the new thin disk, the velocity dispersion ratio a z /<TR 
depends sensitively on the orbital inclination of the infalling satellite. 



could be used to test thick-disk formation scenarios (see also 
iBournaud et^l2009l) . 

Figure[T3]shows the 04 parameter as a function of isophotal 
surface brightness in the V-band for the final thick disks of 
experiment A. This figure shows that the trend of more boxy 
contours at lower surface brightness limits is maintained after 
the growth of the new disk. Interestingly, at a given fi v , the 
contours are slightly less box y in comparison to those pre- 
sented in Figure 8 of Pa per 1 This is likely caused by the 
growth of the new (flatter and more massive) thin disk. Recall 
that Figure [TJ shows that contours at lower fiy are in general 
closer to the center of the system after the induced structural 
contraction. In order to detect the "boxiness" in the exper- 
iments with growing disks one would need to reach surface 
brightness levels beyond /iy > 24-25 mag arcsec -2 , namely 
> 6-7 mag below the central peak of the thick disks after the 
growth of the new disk, /iy ~ 18 mag arcsec" 2 . 

The bottom panel of Figure Q~3] shows for the same exper- 
iment the isophotal surface brightness as a function of the 
isophotal semi-major axis before and after the growth of the 



new disk. The results in this panel suggest that a possible 
cause for the structural change introduced to the vertical sur- 
face brightness profiles in the outskirts of the final thick disks 
(see Figure @|i is the contraction which brings boxy isophotes 
within the region R < 10 kpc. To this end, the presence of 
boxy contours at the outskirts would enhance the separation 
between inner contours (that are more packed together) and 
outer contours (see Figure [T). This would make more no- 
ticeable the difference in the profile shown in Figure |4] (right) 
between \z\ < 1 kpc and |z| > 2 kpc. 

The structure of the initial thick disks is in good agreement 
with observations of external galaxies and the MW in terms 
of thick-to-thin ratios of both scale-lengths and scale-heights 
(see Section 4.2 of lPaper II) . This conclusion was based on the 
assumption that thick-disk structure would not be strongly af- 
fected by the growth of a new thin disk and that this thin disk 
would follow the same distribution as the cold component of 
the thick disk. Overall, the results of the present paper indi- 
cate that thick disks do respond to the growth of an embedded 
thin disk. It is interesting to investigate how these scale ratios 
compare to observations after the formation of the new thin 
disk. In our new experiments, we find the final thick-to-thin 
scale-height ratio to be ~ 0.6/0.125 = 4.8, which is still within 
the range 2.4-5.3 observed in the MW and in SO galaxies (see 
i Buser et all 1 19991: iPohlen et al.1 12004 lYoachim & Dalcantonl 
120061) . 

In what follows we assume that the Galactic thick disk 
originated from the vertical dynami cal heating of a primor- 
dial th in disk by infalling satellites dKazantzidis et al.ll2008b 
iPaper J) . It is interesting to use the results above and attempt 
to determine the initial properties of the initial thick disk of 
the MW, namely those before the structural contraction in- 
duced by the formation of the current thin disk; as well as 
the properties of the pre-existing disk, namely those before 
the interaction with the satellite. For example, we typically 
find that the scale-length of the thick disk after the encounter 
with the satellite increases by a factor /* ert , er ~ 1.3— 1.45, 
as shown in Section 3.4 of IPaper \ On the other hand, the 
formation of a factor of 5 more massive thin-disk compo- 
nent leads to a decrease in the scale-length of f^ om ~ 0.6. 
This implies that the scale-length of the pre-existing disk 
is given by Rf™ = (f« erga . x f? 



) _1 x R d . thlck & or Rf^ 1 



1 Of course, such a relation is only valid under the assumption that the only 
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1 . 1 which is comparable to the present-day value of 
the thick disk. This suggests that the radial expansion caused 
by the accretion event is nearly completely balanced by the 
contraction induced by the growth of the thin disk. More- 
over, using current estimates for the thick-disk scale-length 
of ^d,thick = 2.5 - 3 kpc, we deduce the initial scale-length 
of the Galactic pre-existing disk to be R [ f tiRl ~ 2.9-3.5 kpc. 
Of course, this could represent a considerable problem, since 
disks are expected to hav e been (significantly) smaller in the 
past (Bui trago et"ai]l2008h . 

Similarly, we may compute the initial scale-height of the 
thick disk of the MW from /^ elger ~ 4.7 and f^ om — 0.6 (both 
for a 30° inclination encounter). Thus, after assuming for the 
current thick disk a value of z^thkk = 1 kpc, we compute the 
scale-height of the pre-existing disk as z™ tlal ~ 340 pc. There- 
fore, we estimate the vertical structure of the pre-existing disk 
to be similar to that of the current Galactic thin disk. 

Especially relevant in this context is the study of 
Elmegreen & Elmegreen (2006) who investigated thick disks 
in the Hubble Space Telescope Ultra Deep Field (UDF). 
These authors examined whether the process of adiabatic 
thick-disk contraction could have determined the present-day 
scale height of the thick disk of the MW who se midplane den - 
sity ratio of thick-to-thin disk is 12% (e.g.. lJuric et al.| [2008). 
Their calculations showed that if the present thick-disk com- 
ponent of the MW (with an estimated exponential scale-height 
of 875 pc) began as an equilibrium pure-thick disk at a young 
age, and if subsequent accretion of the entire thin disk was 
adiabatic, then the initial thick-disk sech 2 scale-height had to 
be ~ 3 kpc (i.e. exponential scale-height ~ 1.5 kpc). Such 
a value is considerably larger than that observed for young 
thick disks in the UDF, where the average scale height is 
1.0 ±0.4 kpc. We note that the magnitude of the contrac- 
tion induced in the thick disk by the growth of the new thin 
disk reported by Elmegreen & Elmegreen (2006), f£ om - ~ 0.6, 
agrees well with our estimates. 

4.2. Kinematics 

Paper I showed that when the thick disks are not decom- 
posed into cold and hot components, their mean rotational 
velocities exhibit clear gradients as a function of height (see 
their Figure 20). The magnitude of these vertical gradients 
depends sensitively on the initial satellite orbital inclination, 
with larger gradients corresponding to lower inclinations. 

Figure [14] shows the variations in the mean rotational ve- 
locity of the final thick disks in experiment A as a function of 
galactocentric distance. Interestingly, after the growth of the 
new thin disk the vertical gradients are retained, maintaining 
th eir depe ndence on satellite orbital inclination. As concluded 
in lPaper 3 . the possible detecti on of such vertical gradients in 
the th ick disk of the MW (e.g. lGirard et al.ll2006t llvezic et ail 
2008) would suggest a satellite accretion event with a low or 
intermediate inclination, if the Galactic thick disk was formed 
via the dynamical heating of a preexisting thin disk. Note also 
the progressive decrease of the mean rotation at outer radii for 
larger inclinations. As we illustrate below, this is due to the 
fact that stars from the heated disk and those from the satellite 
have rather different z-components of angular momentum at 
those radii (see Figure [T6i>. 

relevant processes are the satellite accretion event and the thin-disk growth, 
and that secular evolution processes are unimportant in establishing the galac- 
tic structure. 



Paper I also demonstrated that the velocity dispersion ratio 
cf z /(Jr is a fairly accurate discriminator of the initial orbital 
inclination of the accreted satellite (see their Figure 21). Re- 
cent observations in the solar neighborhood report a-. /<jr ~ 
0.6 (e.g. see lLavdenet all Il996t IChiba & Beersl 120011: 
Soubi ran et alJl2003HAlcobe & Cubarsill2005HVallenari et alJ 
2006). Assuming that the Galactic thick disk was formed ac- 
cording to the disk heating scenario, such a value may be sug- 
gestive of an enco unter wi th a satellite on a low/intermediate 
orbital inclination (IPaper J) . Figure Q3] shows the final thick- 
disk dispersion ratio a z j or as a function of galactocentric dis- 
tance in experiment A. It is interesting that the conclusions 
previously drawn are still valid after the growth of the new 
thin disk. An interaction with a satellite on a low/intermediate 
initial inclination is still required to obtain a value of cf z /<jr ~ 
0.6. We note that similar results are valid for the rest of the 
experiments presented in this study. 

For the experiment with initial inclination 30°, we find 
that after the satellite accretion event the thick-disk ve- 
locity dispersions increase by factors (/r, fy, / z ) merg er ~ 
(2,1.7,1.5). These values are obtained by measuring the 
velocity dispersions at 2AR j and at 2.4/^thick ( see Table 1 
and Figure 14 of Paper I, for "z=l"). Additionally, the 
growth of a factor of 5 more massive new disk further in- 
creases the velocity dispersions by factors ( /r, fy,f i)camr ~ 
(1.8,2,2.2), where the final velocity dispersions were mea- 
sured at 2.4/^™j ck (see bottom panel of Figure |5J, As in 
the case of computing the scale-height of the pre-existing 
disk, we can estimate its initial velocity dispersions as 
((T^ m ,af m ,af m ) ~ (18,16,12)kms- 1 . This calcula- 
tion assumes obse rved values for the Galactic thick disk 
(65,54,38) kms" 1 (lLavden et al.lll996t IChiba & Beersll200ll: 
Soubir an et"aT]l2003HAlcobe & Cubarsill2005HVallenari et alJ 
2006; IVeltz et a l.l l2008h . The entire process of heating and 
subsequent contraction leads to an increase by a factor of 
~ 3.4 in all velocity dispersions for this particular example. 

4.3. Phase-Space 

IVillalobos & Helmil (I2009L hereafter [Pajirll showed that 
the distribution of the z-component of the angular momen- 
tum L Z (R) is a good discriminator for separating heated disk 
and satellite stars, especially at large radii (see their Figure 6). 
Figure [16] shows both the initial and final L z distributions of 
disk and satellite stars for experiment A. Stars are located 
within 2 < R < 7 kpc and \z\ < 1 kpc. Interestingly, the 
initial separation between disk and s atellite s tars is now less 
clear compared to that in Figure 6 of IPaper III . This is likely 
due to differences in the intrinsic kinematics of the satellites 
used in each case (spherical versus disky). After the growth 
of the new disk two effects are expected regarding the L z 
distributions of stars. First, disk stars are expected to have 
a steeper slope in their trend of L Z (R), while satellite stars 
should roughly maintain their almost flat slope. This is a con- 
sequence of both the radial contraction experienced by thick 
disks and the adiabatic invariance of L z of each star. Second, 
both type of stars are expected to have a larger L z (= RV^) dis- 
persion at a given radius R. This is because of the increase 
in the V$ velocity dis pers ion measured after the growth of 
the new disk (Section [3.2b . Both effects can be observed in 
the left panels of Figure [16] As the slope of disk stars gets 
steeper, it is expected that the separation between disk and 
satellite stars would become more noticeable. However, the 
significant increase in the velocity dispersion appears as the 
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FIG. 16. — Left panels: Initial and final scatter plots of vertical angular momentum, L z , versus galactocentric distance, R, for disk (grey points) and satellite 
stars (black points). Results are presented for experiment A and for stars that were initially located within 2 < R < 1 kpc and \z\ < 1 kpc. For clarity, only one 
every five stars is plotted. Right panels: Initial and final L- histograms for disk and satellite stars that were initially located within 8 < R < 9 kpc and \z\ < 1 kpc. 
Results are presented for experiment A and dotted lines show the total distributions. 
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FIG. 17. — Histograms of final heliocentric line-of-sight velocities for thick-disk stars after subtracting the mean rotation of the final thick disks. Results are 
presented for experiment A and histograms include disk and satellite stars within a slice around / ~140°, where the contribution of satellite stars to the wings is 
maximal. The width of /-slices is 40°. Dashed lines present the best-fit Gaussian distributions to the histograms. Open points show the fraction of satellite stars 
in each velocity bin, while filled points denote the probability of the observed number of stars compared to what is expected from the best-fit Gaussian in each 
velocity bin. 



dominant effect, making the separation even less evident than 
before the growth of the new disk. 

The right panels of Figure [16] present the distributions of 
disk and satellite stars that were initially located in the solar 
neighborhood (8 < R < 9 kpc) and within \z\ < 1 kpc. The 
histograms show that while there is no clear distinction be- 
tween disk and satellite stars, something that can also be seen 
in the initial thick disks (except in the 60° case), it is possible 
to identify a large contribution of satellite stars in the wings 
of the distributions at low L z . This causes significant asym- 
metries in the total distributions. Nevertheless, it is important 
to note that the trend in L z of disk stars as a function of R 
is maintained and thus should be observable in the Galactic 
thick disk if it was formed via the heating of a pre-existing 
di sk. 

Paper II also showed that the wings of the distributions of 
heliocentric line-of-sight velocities, vi os , measured in small 
volumes in the final thick disks, contain mostly satellite stars. 
These velocity distributions were found to differ significantly 
from Gaussians because of the more heavily populated wings 



(see Figure 10 of iPaper III) . Figure [171 shows a very simi- 
lar behavior for the final thick disks of experiment A. In this 
case, the volumes are located somewhat closer to the center, 
at ~ 3 kpc (1 .5 kpc radius) to take into account the radial con- 
traction of the systems after the growth of the new disk. How- 
ever, we remind the reader that the results were not found to 
depend on either the size or the location of the volumes. It is 
also important to stress that the dispersions of the Vi os distribu- 
tio ns are sig nificantly larger than those presented i n Figure 10 
of Pap er 11 Nonetheless, the same tests used in Pap er III to 
quantify the statistical significance of features present in the 
vios distributions can successfully detect the contribution of 
satellite stars to the wings in the present study. These statis- 
tical tests measure the likelihood of finding satellite stars at a 
given velocity bin of a vi os distribution by generating a num- 
ber of random realizations based on the best Gaussian fit to 
the distribution. 

5. SUMMARY AND CONCLUSIONS 

Using a suite of collisionless A^-body simulations we have 
examined the evolution of the structural and kinematical prop- 
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erties of simulated thick disks induced by the growth of an 
embedded thin disk. Our simulation campaign quantifies the 
importance of various parameters related to the growing disks 
that could influence the response of thick disks, including: (1) 
the growth timescale; (2) the final mass; (3) the vertical scale- 
height; (4) the radial scale-length; and (5) the sense of rota- 
tion of the initial thick disk with respect to its halo (prograde 
or retrograde). 

Our simulations demonstrate that after the growth of the 
new thin disk, thick-disks experience a strong radial as well 
as vertical contraction which leads to a significant decrease in 
their scale-lengths and scale-heights. This contraction is trig- 
gered by the deepening of the potential well of the system due 
to the mass growth associated with the new disk. An impor- 
tant consequence of this contraction is related to the migration 
of thick-disk stars from the outskirts inwards. Stars that were 
typically at R ~ 15 kpc before the growth of the new disk, are 
found at R ~ 9 kpc after the growth is completed. 

Despite the strong contraction, the mass fraction associated 
with the kinematically cold component of the final thick disks 
remains < 25%. This is comparable to the mass prese nt in the 
remnants of the satellite accretion events discussed in lPaper A 
The outskirts of these disks (R ~ 3.8-5.6^ .thick) present boxy 
contours at very low surface brightness levels (> 6-7 mag 
below the central peak). 

We find that the thick-disk structural and kinematical evo- 
lution is driven primarily by the total mass and scale-length of 
the growing thin disk, with a weak dependence on the char- 
acteristics of the initial thick disk. Conversely, the thin-disk 
growth timescale appears to have a minor influence on the 
final thick-disk structure provided that this timescale is suffi- 
ciently long for the process to be considered adiabatic. Sim- 
ilarly, neither the thinness (in terms of scale-height) nor the 
alignment of its angular momentum vector (either to the halo 
angular momentum or to that of the thick disk) seem to affect 
the final properties of the thick disks. The sense of rotation of 
the initial thick disk has an effect on the radial extension of the 
cold component of the final thick disk, which is significantly 
smaller for retrograde rotation. 

Kinematically, the growth of a more massive thin disk in- 
creases the mean rotation and the velocity ellipsoid of the ini- 
tial thick disks. The increase in the mean rotation is associated 
with the larger mass deposited at a given radius after the struc- 
tural contraction. Moreover, as a consequence of the contrac- 
tion and additional mass accumulation on the midplane of the 



thick disk, the velocity dispersion increase is linked to a sig- 
nificant steepening of the gravitational potential gradient near 
the disk plane. 

We close with a few words of caution and a discussion 
of fruitful directions for future work that may lead to more 
conclusive statements about the role of a reforming thin disk 
in establishing the structural properties of thick disks. We 
reiterate that we have only modeled the gravitational po- 
tential of a thin disk that slowly grows in mass over time 
in the collisionless regime, and we have neglected both the 
cosmological context and the complex interplay of effects 
(e.g., gas cooling, star formation) relevant to the formation 
and evolution of spiral galaxies. Previous studies have been 
carried o ut in the context of galact i c mergers including gas 
physics (iBarnes & Hernauisl Il996t iKazantzidis et al] 120051: 
I di Matteo et all 120071: iHopkins et alJ l2008t iBournaud et alJ 
120091: fCallegari et alj|2009t iMoster et al Il2009t) and within a 
cosmological context dGovernato et alJ 120071) . therefore we 
consider that an important future study will be to follow self- 
consistently, as the satellites fall in, the growth of the new thin 
disk via accretion-induced bursts of star formation, deposition 
of gas by the accreting satellites themselves, and smooth gas 
accretion through the cooling of hot gas in the galactic halo. 
While a full exploration of these contingencies is challenging, 
we intend to extend the present investigation in this direction 
in a forthcoming work. 
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